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Recent progress in the simulation of Casimir forces between various objects has allowed tradi- 
tional computational electromagnetic solvers to be used to find Casimir forces in arbitrary three- 
dimensional objects. The underlying theory to these approaches requires knowledge and manipula- 
tion of quantum field theory and statistical physics. We present a calculation of the Casimir force 
using the method of moments via the argument principle. This simplified derivation allows greater 
freedom in the moment matrix where the argument principle can be used to calculate Casimir 
forces for arbitrary geometries and materials with the use of various computational electromagnetic 
techniques. 



INTRODUCTION 



DERIVATION 



The ability to model the bulk effect of molecular forces 
has become increasingly important with the development 
of micro- and nanoelectromechanical systems (MEMS 
and NEMS). The static and dynamic behaviors of such 
devices involve the evaluation of molecular and electrody- 
namic forces. The Casimir force can simulate them more 
accurately than a simpler form of the van der Waals force. 
As such, a simple and robust method of calculating the 
Casimir force is desirable. Recent methods have been de- 
rived based around boundary element methods PQI2]. The 
Reid, Rodriguez, White, and Johnson (RRWJ) method 
allows the use of computational electromagnetic (CEM) 
solvers using the method of moments to find the Casimir 
force for arbitrary three-dimensional objects of arbitrary 
medium. This presents an attractive method for cal- 
culating the Casimir force. However, this method uses 
a traditional method of moments matrix which suffers 
from low frequency breakdown and mesh density break- 
down. Instead of deriving through the path integral, we 
show a new method that uses the argument principle with 
boundary element methods to arrive at an improved for- 
mulation. 

Calculating the Casimir force using the argument prin- 
ciple was originally applied by Van Kampen et aZ.[5] for 
a much simpler problem. With the new derivation, we 
can choose from a wide range of modern CEM formu- 
lations that overcome shortcomings with the minimum 
of work. In doing so, we expand beyond the class of 
problems that can currently be solved using the RRWJ 
method. A simple demonstration of this is done using 
a formulation that avoids the low frequency breakdown 
for perfect electrical conductors (PEC), the Augmented- 
Electric Field Integral Equation (A-EFIE) 4 approach. 



The RRWJ method starts out using a path-integral ex- 
pression for the Casimir energy of arbitrary PEC objects 

at r = o[H|S]. 
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Here, the functional integration Z is performed over all 
possible configurations of the current J(r) on the surfaces 
of the objects and G K (r,r') is the dyadic Green's func- 
tion with the Wick rotated imaginary frequency. Z^ is 
computed similarly as Z except with all objects removed 
to infinite separation. By expanding the continuous cur- 
rent distribution J(r) over a finite discrete set of basis 
functions f m , defined over the surfaces of the objects, the 
RRWJ method reduces to 
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The elements of the matrix M, being 



f i (r).G K (r,r , ).f i (r , )dr'dr (5) 



are of the same form as the matrix elements of the tra- 
ditional CEM impedance matrix for the method of mo- 
ments. Equation Q is derived from Equation ([I]) assum- 
ing a specific action for the path integral and properties 
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of the matrix M. Attempts to modify the resulting im- 
pedance matrix used in Equation Q requires a new path 
integral formulation as well as the assumption that the 
matrix M be a positive or negative definite matrix. In 
addition to PEC objects, this formulation can handle ho- 
mogeneous dielectric objects[5]. 

However, greater freedom is needed to have a matrix 
M that handles a wider class of problems, like inhomo- 
geneous dielectric objects, or addresses numerical insta- 
bilities like low frequency integration points. To ensure 
low frequency stability, the loop-tree decomposition^ 
can be used, whereby the impedance matrix is the re- 
sult of a similarity transform of the traditional method 
of moments impedance matrix. This transform, achieved 
using left and right matrix multiplications, is canceled 
out when the determinant is normalized with Moo, thus 
resulting in the mathematical equivalence of Equation 
Q. Other low frequency stable methods like A-EFIE, 
unfortunately, are not symmetric positive definite nor 
represent a similarity transform of the predefined ma- 
trix M. Instead of evaluating a new path integral, we 
approach the problem using the argument principle. 

Note that the Casimir force results from the relative 
perturbation of the quantum vacuum fields by the perti- 
nent objects. The objects force the fluctuating vacuum 
fields to conform to the appropriate boundary conditions 
and, in doing so, change the energy density of the vac- 
uum. If one determines the eigenfrcquencies of the field 
configurations that satisfy the geometry's boundary con- 
ditions, then the unnormalized energy of the geometry 
can be found via [71 |S] 
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Previous derivations of the Casimir force have made use 
of the argument principle to find the eigenmodes of the 
vacuum fields to define the energy density [51 [71 |9"HTT] . To 
calculate the Casimir energy from these dispersion equa- 
tions, we first note that the argument principle states 
that 
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where loq^ are the zeros and Wooj are the poles of the 
function f(ui) inside the contour of integration. Noting 
that the Casimir energy of a geometry is given by Equa- 
tion ([6]), we can relate the Casimir energy with the con- 
tour integral from above. Using integration by parts, the 
Wick rotation where u — ick, and folding the integral, 
we find 
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where the zeros of /(w) correspond to the eigenmodes of 
the structure. It is necessary, in calculating the Casimir 



energy, to renormalize the vacuum energy by subtracting 
off an appropriate normalization energy represented by 
modes uj n . 



That is, 
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Following the same process, the Casimir energy is repre- 
sented as 
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The theoretical problem is then to find some function 
f(uj) and its associated normalization /n 0rm (w), taken 
as the geometry when the objects are infinitely sepa- 
rated, that are zero at the eigenfrequencies of the Casimir 
problem. This has been previously done by finding the 
dispersion relationship of the geometry using primarily 
closed-form analysis. This could only be applied to a 
small number of geometries where exact dispersion rela- 
tions could be derived like periodic gratings or dielectric 
slabs and was not extended to arbitrary geometries. 

The preceeding derivation relies on the physical inter- 
pretation of the natural modes of the problem and the 
summation over their energies to arrive at the Casimir 
energy in Equation As such, the above is valid for a 
lossless system where the eigenfrequencies are real, but 
the physical meaning of these eigenfrequencies becomes 
indeterminate in a lossy system with complex eigenfre- 
quencies. Despite this, the use of the argument principle 
using lossy materials still seemingly resulted in the cor- 
rect Casimir force [SI HH H3J. Only recently have papers 
have been published that rigorously explain whv |14lHB] . 

The traditional CEM formulation for PEC objects is 
derived assuming that the fields arising from currents in- 
duced on the surface of a PEC must cancel any excitation 
fields that are present. This boundary condition leads to 
the Electric Field Integral Equation (EFIE), 

-h x E s (r) = n x E l (r) = hxiujfi J dr'G (r, r') • J(r') 

(11) 

where G (r, r') is the non-Wick rotated dyadic Green's 
function, E s is the scattered electric field and E l is the 
incident exciting electric field. The method of moments 
is used to derive an impedance matrix to relate surface 
currents to the excitation fields. That is, 



Z J = V 
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where Z is the matrix representation of the real frequency 
Green's operator and J and V are the vector representa- 
tions of the current J(r) and exciting fields respectively. 
As previously stated, the eigenmodes of the quantum vac- 
uum are those where the fluctuating fields satisfy the 
boundary conditions of the geometry. In other words, 
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eigenmodes exist where we have currents on the surface 
of the objects without any excitation fields. These natu- 
ral modes automatically satisfy the geometry's boundary 
conditions, and they also satisfy the relation, 



Z J = 



(13) 



Being a singular matrix, we conclude that / = det Z = 
for the natural eigenfrequencies where the vacuum 
fields satisfy the boundary conditions without an exter- 
nal source. Artificially enclosing the objects in a PEC 
cavity ensures that the eigenmodes of the system will 
be lossless and discrete. By using an appropriate cavity 
dyadic Green's function in Z, dct Z can be used as the 
dispersion relation in Equation (Jsj) to find the energy of 
the bounded cavity. Expanding the volume of the cav- 
ity to infinity and applying the normalization where the 



objects are infinitely separated, we reach Equation (10). 

However, by using the cavity Green's function, the re- 
sult is still incorrect and thus an infinitesimal loss, to 
enforce the Sommerfeld radiation condition [17], is added 
to remove the reflections off the cavity walls at infin- 
ity. In doing so, the resulting Green's function becomes 
the free-space Green's function, as it was in the origi- 
nal EFIE, and the RRWJ formulation is rederived. Due 
to the infinitesimal loss, the eigenfrequencies of the EFIE 
are complex frequencies requiring one to justify its deriva- 
tion for lossy systems. The strength of this derivation is 
its large degree of flexibility. As long as one finds a rela- 
tionship that similarly relates the geometry's boundary 
conditions with the fields excited by induced and external 
sources, then the same procedure can be applied. This al- 
lows us to rederive the RRWJ formulations for PEC and 
dielectric homogeneous objects in addition to more ad- 
vanced CEM formulations that would otherwise require 
novel evaluations of the original path integral. 

One difficulty that arises with the numeric integration 
of the Casimir energy and force over the wavenumbers 
k is the low frequency breakdown. Looking at the force 
spectrum for two PEC rounded cylinders, as shown in 
Figure [T] the spectrum is concentrated in the lower fre- 
quencies as is expected due to the normalization. Evalu- 
ating the eigenvalues for the impedance matrices at lower 
frequencies becomes inaccurate due to the low frequency 
breakdown of the EFIE impedance matrix. While in- 
creasing the precision of the calculations helps address 
this issue, it is still a problem depending upon the CEM 
formulation and class of problems being solved. A bet- 
ter remedy is to use a more appropriate formulation for 
the impedance matrix, like A-EFIE, that will not suf- 
fer from the same breakdown. A-EFIE starts from the 
normal EFIE formulation whose matrix representation in 



Equation ( 12 ) can be equivalently expressed as 
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Assuming that the subspace is spanned by the Rao- 
Wilton-Glisson (RWG) basis Q3], the problem has e RWG 
edges and p patches. The vector and scalar potential 
matrices above, V and S respectively, and an additional 
matrix representation of the scalar Green's function us- 
ing patch basis are defined to be 



V r , 
S, 



/A 



/ A m (r)- f g(r,r')A n (r')dS'dS (15) 



s„ 



, , 1 / V ■ A m (r) • / g{r, r')V • A n {r')dS'dS 

(16) 

Mr) / 5 (r,r')Mr')dS'dS (17) 



where g(r, r') is the homogeneous scalar Green's function, 
A(r) is an RWG basis that is not normalized against the 
edge length, and h(r) is a pulse basis. We can relate 
the patch matrix, P, with the scalar potential matrix, S, 
using the incidence matrix D by 



S = D P D 

By accounting for current continuity, 

D • J = ikgcop 
we can derive the A-EFIE matrix system as 
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The given formulations for A-EFIE can be further im- 
proved upon and it is seen that the A-EFIE impedance 
matrix, Z^, will have the function f(u) = det Z A = at 
the eigenfrequencies of the geometry. The A-EFIE im- 
pedance matrix can be used as a direct replacement for 
the EFIE impedance matrix in the Casimir calculations 
by use of the argument principle. Thus, the Casimir en- 
ergy and force become, 
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where Vj represents the derivatives with respect to a 
physical displacement of the z-th object (the one we wish 
to find the forces acting upon). The integrand of the 
Casimir energy is found by solving for the eigenvalues of 
Z A and Z Aoo . 
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where A„ and A£° are the eigenvalues of Za and Za,<x> re- 
spectively. For the Casimir force, note that the integrand 
can be expressed as 
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where a n are the eigenvalues of the generalized eigenvalue 
problem 



V^Z^ • x = <xZa ■ x 



(25) 



It should be noted that the gradient acting upon Za acts 
on all elements of the matrix, making the bottom block 
matrices zeros. That is, 



first time, the low frequency A-EFIE was used to over- 
come the limitations of the EFIE formulation in Casimir 
force calculations. This opens up possibilities for other 
improvements like using the mixed-form fast multipole 
algorithm with A-EFIE [H Q15] or the use of other for- 
mulations like the PMCHWT method for homogeneous 
dielectrics [2"Dl - f2"5] (which can be used to rederive another 
result by RRWJ[2]), the volume integral method for inho- 
mogeneous dielectric objectsJU], and the layered medium 
Green's function|25j to help model systems with infinite 
slabs like a substrate. Many of these formulations open 
a new class of geometries or scale of problems that would 
otherwise not be accessible using the boundary element 
methods found in current methods. 
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Using A-EFIE, we can more accurately calculate the 
low frequency energy spectrum than using EFIE while at 
the same time being able to solve the same set of prob- 
lems as the original RRWJ formulation. While in gen- 
eral, double precision can give satisfactory results, Figure 
[Q shows that the A-EFIE and EFIE results start to dif- 
fer as the frequency drops when using single precision for 
the force between two parallel PEC rounded cylinders of 
length to radius ratio of 6, and separation to radius ra- 
tio of 4. The nominal dimensions are immaterial as the 
force and energy results are normalized with respect to 
the radius of the objects. The result from the numer- 
ical integral is calculated using points indicated in the 
figure and the points with the largest contribution lie in 
the region of poor performance. The deviation in the 
calculated integration is a result of the location of the 
Gaussian points but the error in the EFIE result can be 
greater than 100% while A-EFIE remains in agreement. 
Figure [2] demonstrates how the total force between two 
PEC spheres can diverge over a range of values for the 
ratio between the separation and radius of the spheres. 



CONCLUSION 

This paper, using the argument principle, presents an 
alternative and novel derivation for the Casimir energy 
and force using boundary element methods. As a re- 
sult, a simple procedure can easily be adapted to dif- 
ferent CEM methods to derive previous formulations or 
new ones. For example, the characteristics of the Casimir 
force can present difficulties when using traditional EFIE. 
The importance of the low frequency content of the in- 
tegrand gives rise to erroneous results when relaxing the 
computational precision due to low frequency breakdown. 
As a demonstration of the utility in using boundary 
element methods with the argument principle, for the 
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FIG. 1. Comparison of the integrand for the Casimir force 
integral between capsules for the EFIE and A-EFIE formula- 
tions using single precision. 
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FIG. 2. Comparison of the Casimir force between two PEC 
spheres for the EFIE and A-EFIE formulations using single 
precision. 



